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ABSTRACT 


The population of Low Surface Brightness (LSB) galaxies, which are objects 
with central surface brightnesses at least one magnitude fainter than the night 
sky, is crucial for understanding the extremes of galaxy formation and evolution 
of the universe. As LSB galaxies are mostly rich in gas (Hi), the a.dO-SDSS DR7 
sample is absolutely one of the best survey combinations to select a sample of 
them in the local Universe. Since the sky backgrounds are systematically over¬ 
estimated for galaxy images by the SDSS photometric pipeline, particularly for 
luminous galaxies or galaxies with extended low surface brightness outskirts, in 
this paper, we above all estimated the sky backgrounds of SDSS images accu¬ 
rately in both g and r bands for each galaxy in the a.dO-SDSS DR7 sample, using 
a precise method of sky subtraction. Once subtracting the sky background, we 
did surface photometry with the Kron elliptical aperture by the SExtractor soft¬ 
ware and htted geometric parameters with an exponential prohle model by the 
Galht software for each galaxy image in the a.dO-SDSS DR7 sample. Basing on 
the photometric and geometric results, we further calculated the B-band central 
surface brightness, /io(B), for each galaxy and ultimately dehned a sample of LSB 
galaxies consisting of 1129 galaxies with po(B) > 22.5 mag arcsec”^ and the axis 
ratio b/a > 0.3 from the 12423 a.dO-SDSS DR7 galaxies. This Hl-selected sample 
of LSB galaxies is a relatively unbiased sample of gas-rich and disk-dominated 
LSB galaxies, which is complete both in Hi observation and the optical magni¬ 
tude within the limit of SDSS DR7 photometric survey. This LSB galaxy sample 
spans from 22.5 to 28.3 in /io(B) with a fraction of 4% fainter than 25.0 mag 
arcsec”^ in B-band central surface brightness, and distributes from -27.0 to -12.3 
mag in the absolute magnitude in B band (M(B)), including 43 faintest galaxies 
(M(B) > -17.3 mag). This sample is an absolutely blue LSB galaxy sample, of 



which 98% are bluer than 0.75 mag in B-V color. As for 21cm Hi properties, our 
LSB galaxy sample has a median =0.87, and a large portion (95%) has 

high mass of Hi (M(Hl) > 10^'^M©), which supports that galaxies in this LSB 
galaxy sample are mostly gas-rich. Additionally, we statistically investigated the 
environment of our LSB galaxies, and found that up to 92% of the total LSB 
galaxies have less than 8 neighbouring galaxies, which strongly evidenced that 
LSB galaxies prefer to reside in the low-density environment. 

Subject headings: catalogs - galaxies:spiral - galaxies:irregular -techniques:photometric 
-methods:data analysis - methods:statistical 
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Introduction 


Among hundreds of billions of galaxies in the universe, there exists such a class of 


galaxies called low surface brightness (LSB) gala xy which is at leas 


than the night sky in central surface brightness fjimpey fc Rothun 


one magnitude faint er 


19971: 


Vollmer et ah 


2 OI 3 I ). In most cases, the central region of a galaxy is the brightest part of the whole 


galaxy. Thus, a galaxy with its B-band central surface brightness fainter than a certain 
threshold value is traditionally regarded as a LSB galaxy. Howeyer, this threshold value is 


not unihed up to now but it usu a^ 


papers (e.g. 


Impev et ah 


feooil); 


ly varies between 22.0 and 23.0 mag arcsec ^ in different 


Ceccarelli et ah 


(120121) 1 . Besides B band, central surface 


brightnesses in optical red or even near-infrared band have also been adopted to attempt 
to distinguish LSB galaxies from high surface brightness (HSB) g alaxies, for exarnple. 


R-band central surf ace brightness fainter than 20.8 mag arcsec ^ (iCourteau et ah 


1996 


Adami et ah 


2006h. or Ks-band central surface brightness fainter than 18.0 mag arcsec 


-2 


flMonnier Ragaigne et ah 1120031 ) 


Limited by the current capacity of observation, LSB galaxies can not be quite easily 
observed from ground-based optical telescopes due to the faint nature of LSB galaxies 
themselves and the comparable sky background. Therefore, contributions of LSB galaxies to 
the universe have been underestimated for a long time, which may further lead a bias to our 
understanding of the universe. Due to missing of galaxies with faint surface brightness as 
a result of the surface brightness selection effect in observation, the faint end of the galax y 


luminosity 


McGaugh 


funct ion should be seriously underestimated before flDisney &: Phillipps 


1983 


19961 ). Therefore, LSB galaxies would pl ay an important role in improving the 


study of the faint end of galaxy luminosity function (jSprayberrv et ah 


1997 


Blanton et ah 


2005 


Geller et ah 


2 OI 2 I) . Besides to the faint end study of galaxy luminosity function, LSB 


galaxies are proposed to be significant in many other aspects. For example, they are helpful 
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for us to study the star formation process in the low gas surface dens ity flNoguchi 


the baryonic matter density and the galaxy formation and evolution fjO’Neil et al 


2001 


20001 ), 


and also usef ul for checking 


the Universe (jPeebles et ah 


he pr ediction of the dark matter to the large-scale struction of 


200 ih . 


Recently, with the modern techique of observation rapidly improved, more LSB galaxies 
in the Universe, particularly in the local universe, can be detected. More works have started 
to focus on studies of the LSB galaxies and found the signihcance of this type of galaxies on 
improving our understanding of the Universe. For instan ce, SPSS is on e of th e databases 


well suited for searches for and studies of LSB galaxies. 


Kniazev et al. 


( 20041) identihed 


129 LSB galaxies with /io(g) down to 25.0 mag arcsec ^ from SDSS Early Data Release 
(EDR) using their own developed method, among which several unusal LSB galaxies were 


found, such as one with a structure t 


LSB galaxy candidates, 
disk LSB galaxy (/io(B) 


r at loo ks like two exponential disks and ten new giant 


Liang et al. fl2007l ) serendipitously found a new nearby edge-on 
~ 23.68 mag arcsec“^) with low metallicity (12 -|- log(0/H)~ 


7.88) from the SDSS, which was mistakenly classihed as a star in the SDSS DR4 database 
because the “star” was actually one of its Hll regions hosted by galaxies with intrinsically 
low global luminosities and low surfa c e brig htnesses. Then, from the main galaxy sample 


of SDSS DR4 database 


Zhong et al. 


(120081 ) established a large sample of 12282 nearly 
face-on LSB galaxies with /ro(B) between 22.0 and 23.6 mag arcsec“^. Based on SDSS data 
of this LSB sample, spectroscopic properties including dust extinction, strong emission-line 
ratios, oxygen abundan ces, nitrogen- t o-oxy gen abundance ratios of LSB galaxies have been 


statistically studied by 


Liang et al. 


(120101) . which showed that LSB galaxies with lower 


surface brightness generally have lower metallicities, dust extinction and stellar mass. 
Compared with high surface brightness (HSB) galaxies, it was suggested that LSB galaxies 


have not had dram atically different s 


;ar formation and chemical enrichment histories 


from HSB galaxies (iGao et al. 


2010 


Liang et al. 


20101) : however, the red and blue LSB 
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galaxies perhaps have different star formation histories:blne LSB galaxies are more likely 
to experience a sporadic star formation event now while r ed LSB galaxies are more likely 


to ha ve formed stars continnously over the past 1-2 Gyr flGao et al 


ml 


Zhong et ah 


2 OIOI ). It is found that LSB galaxies tend to have a lack of companions compared to HSB 


galaxies at small scales and the isolation is more connected with the survival, formation 
and evolution of LSB galaxies. When compared to a sample of isolated LSB galaxies, LSB 
galaxies which are close to neighbors clearly show relati vely high star form ation signatures 


and/or have a high population of recently formed stars flGalaz et ah 


2 OI 1 I). 


Recent studies indic ated that LSB galax ies contribute 20% to dynamical mass of 


galaxies i n the universe ( 


2006 


galaxies flMcGaugh 


Vlinchin et ah 


1996 


Haberzettl et ah 


Bothun et ah 


2004 ). 30% ~ 60% to the number density of loc al 


1997 


O’Neil et ah 


2000 


Trachternach et ah 


20071 ). which suggests that the contribution of LSB galaxies to the 


universe, especially to the local universe, can not be negligable. Although, Malin-1, the 
hrst discovered LSB galaxy, is a giant spiral galaxy, LSB galaxies are in fact dominated 
by faint galaxies because as high as 95% of the mass of LSB galaxies are occupied by 
non-luminous dark matter. They are morphologically extended, disk-like or irregular, and 
mostly rich in gas, poor in dust, and deficient in metal, and the chemic al evolution of LS B 


galax ies differs from that of galaxies that define the Hubble Sequence (jPustilnik et al. 


201 ih. The oxygen abundances of LSB galaxies reac 


ti as 


flMcGaugh fc Bothun 


1994 : 


Roennback fc Bergvall 


1995 : 


ow as 1/3 1/5 so^ 


Burkholder et al. 


or abundance 


2011 ). and so 


far, the already known galaxies with the lowest metallicities belong to LSB galaxy type. 


Proved bv practical observations. LSB galaxies have a very smal 

number of Hll regions. 

whic 

1 are also small in size and weak in emission ( 

Schombert et al. 

1992; 

McGaugh et al. 


1995 


Schombert et al. 


20131 ). This observational phenomen on indicates that LSB 


galaxies have very low star formation rates (SFRs) at present (Ivan der Hulst et al. 


1993 
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van den Hoek et al. II 2 OOOI: iGalaz et al. 1 1201 ill) while they are still abundant in gas material 


of star formation. Obviously, LSB galaxies are unevolved and the star fo rmation ac 


of them have be en suppressed at a very low level especially in the disk flMo et ah 


Das et ah 


ivities 


1998 


20091) ■ Although LSB galaxies have low levels of current star foramtion activities, 
they contain variaties of stellar populations from very red (B-V ~ 1.2 mag) to very blu e 


(B-V ~ 0.2 mag) in co lor and nearly covers all the H-R diagram fjO’Neil et ah 


Zackrisson et ah 


1997bl : 


20051) . Up to now, the reasons why LSB galaxies have low star formation 


activities in the condition of still being rich in gas material are still an open question. 

LSB galaxies are important populations in the universe, so we nee d more good s a mples 


to better study them. Although the LSB galaxy sample established by 


Zhong et al. 


(120081) 


mentioned above has greatly extended the known sample of LSB galaxies at that time, 
galaxies with B-band absolute magnitude greater than -18.0 mag were cut out from the 


SPSS DR4 database before LSB galaxy selection. Therefore, the sample of 


Zhong et al. 


(120081) is mainly composed of intermediate LSB galaxies and inevitably lacks the dwarf LSB 
ones in the SDSS database. As LSB galaxies are mostly gas-rich, surveys of gas-rich galaxies 
would provide us an opportunity to select a large number of LSB galaxies, which avoids 
the surfa ce brightness ef f ects th at limit the usefulness of optical surveys for finding LSB 


galaxies. 


Minchin et al. 


(120041) assembled a Hl-selected sample of 129 LSB galaxies from a 


very deep survey for neutral hydrogen (HIDEEP) with the Parkes multibeam system, and 
estimated the cosmological signihcance of LSB galaxies which provide more than half of the 
gas-rich galaxies by number, about 30 per cent of the contribution of gas-rich galaxies to the 
Hi density of the Universe, about 7 per cent of the contribution to the luminosity density 
of the Universe, about 9 per cent of the contribution to the baryonic mass density of the 
Universe, and about 20 per cent of the contribution to the dynamical mass density of the 
Uni verse. The ALFALFA surve y (Giovanelli 2007) is one of the successful 21cm Hi surveys. 


and 


Trachternach et al. 


(120061) performed a blind optical follow-up observations of a part of 





































the region covered by the blind Arecibo Hi Strip Snrvey (AHISS 


Zwaan et ah 


fll997h h and 


detected optical connterparts of all Hi detections, of which 30 per cent were LSB galaxies 
at the limiting snrface brightness of /iB^ijm=25.2±0.31 mag arcs ec"^. Fortnnately, th e data 


of the 40 per cent sky area of the total ALFALFA snrvey (a.40 


Havnes et ah 


(201l|)) have 


been released to the pnblic, and the overlap between a.40 and the SDSS (a.40-SDSS) would 
be undoubtfully one of the best survey combinations for us to search a relatively unbiased 
sample of LSB galaxies. In this paper, we will assemble an order-of-magnitude larger sample 
of LSB galaxies from the a.40-SDSS data, this Hl-selected LSB galaxy sample will be not 
only complete in Hi observation but also in optical magnitude within the observational limit 
of SDSS photometric survey, and it would definitely provide a good (relatively unbiased) 
sample for us to better study the LSB galaxies aligning the SDSS optical, ALFALFA Hi 
bands and other bands. 


As the hrst paper of our series of papers on LSB work, we will search a Hl-selected 
sample of LSB galaxies from the overlap between a.40 and the SDSS DR7, and make 
statistical studies about the physical properties of LSB galaxies of our sample, including 
the optical, radio and environmental properties. Then, for the next few papers we will 
base on this LSB galaxy sample and study the oxygen abundance, the nitrogen-to-oxygen 
abundance ratio and perhaps further discuss the origin of the nitrogen element. Combining 
with the GALEX UV data, the hve-color SDSS optical data and the WISE infrared data, 
and the Ha images of the LSB galaxies of our sample, star formation rates of the LSB 
galaxies will be estimated and we will further study stellar populations to reveal their star 
formation histories and evolution which is one of the important and still open questions 
about LSB galaxies. 

In this paper, we will give a brief introduction of the ALFALFA survey, a.40 catalog, 
SDSS DR7, and the a.40-SDSS DR7 sample in Section 2. The data reduction processes 
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done by ourselves for each galaxy in the a.dO-SDSS DR7 sample will be described in detail 
in Section 3, including sky subtraction (Section 3.1), surface photometry by SExtractor 
(Section 3.2), geometry by Galht (Section 3.3), and calculation of central surface brightness 
(Section 3.4). In Section 4.1, we select out a relatively unbiased sample of LSB galaxies 
from the a.40-SDSS DR7 parent sample according to the selection criteria. Then in Section 
4, we statistically make studies on optical properties (Section 4.2.1), Hi properties (Section 
4.2.2) and the environmental properties (Section 4.2.3) for our LSB galaxy sample. Finally, 
we summarize our work in Section 5. 


DATA 


2.1. The a. 40 catalog 


The Arecibo L-band Feed Array (ALFA), a L-band (1.4GHz) receiver in Arecibo, 
the world’s most sensitive radio telescope at Arecibo Observatory in Puerto Rico, enables 


such an extragalactic H i survey as A 


Giovanelli et ah 


(I2nn5ah : 


Giovanelli 


FALFA (the Arecibo Legacy Fast ALFA survey 


(I2nn7h : 


Haynes 


(120071 ) . Exploiting the large collecting 


area of the Arecibo antenna and its relatively small beam size (~3.5'), ALFALFA is a very 
wide area (7000 deg^ of the sky at high Galactic latitudes) blind extragalacitic Hi survey, 
which aims to catalog all gas—bearing extragalactic objects in the local universe and has 
conducted a deep and precise census of the local Hi universe over a cosmologically signihcant 
volumn. Initiated in February 2005 and completed in 2012, ALFALFA has obtained a Hi 
line spectral database of more than 30000 extragalactic Hi line sources c overing the redshi ft 


range between -1600 and 18000 km s ^ with a resolution of ~5 km s^ flGiovanelli et ah 
2005bll . 


Now a catalog of Hi detections covering about 40% of the full ALFALFA survey sky 
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area, has been publicly released as the a.40 catalog flHavnes et ah Il201l[) . The sky areas 


contained in the a.40 catalog are regions 07^30™ < R.A. < 16^30™', +04° < deck < +16°, 
and +24° < deck < +28° (the “spring” region) and 22^00*" < R.A. < 03^00™, +14° < deck 
< +16°, and +24° < deck < +32° (the “fall” region). The a.40 catalog consists of 15855 
Hi detections, 15041 of which are certainly extragalactic objects and the left 814 are more 


likely to be Galactic h igh velocity cloud 


reservoirs of gas (Hll) fjSchombert et ah 


HVC ). LSB galaxies are mostly containing large 


20131) . so the a.40 catalog provides us an excellent 


database to hunt a large number of LSB galaxies in the local Universe. 


2.2. The SDSS DR7 


Mapping one 


Gunn et ak 


Strauss et ak 


(1998h 


q uarter of the entire sky, the Sloan Digita l 


York et ak 


fl2000h : 


Lupton et ak 


(I2Q01I); 


Sky Survey (SDS S: 


Stoughton et ak 


fl2002h : 


(1200211 ) aims to obtain GGD imaging in hve broad bands (u,g,r,i,and z) and 


spectroscopy from 3800 to 9200Aof mil 


wide-held 2.5m telesco pe (IGunn et ak 


Data Release 7 (DR7; 


ions o f galaxies, quasars and stars, using a dedicated 


Abazaiian et ak 


20061) a t Apache Point Observatory in New Mexico. 


(120091) 1 is the seventh major data release and 


provides images, spectra, and scientihc catalogs. It is the hnal data release of SDSS-II, an 
extension of the original SDSS, which share some footprints with the a.40 data set. 


2.3. The Parent sample:a.40-SDSS DR7 sample 

By cross-referencing the a.40 and SDSS DR7 photometric data sets where the two share 
footprints, the ALFALFA team provides the cross-identihcations of a.40 Hi sources with the 
photometric and spectroscopic catalogs associated with the SDSS DR7 in the a.40 catalog. 
Of the total 15041 extragalactic objects from the a.40 catalog, there are 12468 sources 













































11 


havin g optical counterparts (OCs) in the SDSS DR7 photometric data set flHavnes et al 
20111) -Here, we have to introduce 2 SDSS photometric catalogs, the PhotoObjAll and the 
PhotoPrimary catalogs. The PhotoObjAll catalog contains all measured parameters for all 
photometric objects of SDSS. As an object may be observed two or more times in SDSS 
imaging due to the overlaps at many levels of the imaging (runs, stripes), the PhotoPrimary 
catalog has been created to only contain the best observation of an object. Usually, we 
should take objects in the PhotoPrimary catalog for convincing scientihc studies. Among 
the 12468 ALFALFA sources with SDSS DR7 OCs, there are 12423 sources belonging to 
the PhotoPrimary catalog. So for a convincing scientihc study in this paper, we will only 
regard these 12423 PhotoPrimary sources as our parent sample. 

Here we just give a brief description of this parent sample. As shown in Figure [H the 
parent sample covers a distance range from 0 to 260 Mpc with 84% between 50 and 220 
Mpc, a Hi mass ra nge from 10^’^^ to 10 ^°-^^ Mq with only 1.6% as low Hi mass sources 


(M(Hl)< 10^’^ Mq; iHuang et al. fl2012l )). a heliocentric radial velocity range from -400 
to 18000 km s“^ with 75.7% between 5000 and 15000 km s“^, a Hi line width range from 
9 to 885 km s“^ with 99.7% below 600 km s“^, a r-band magnitude range from 9.73 to 
21.88 mag with 4.5% fainter than 18 mag, and a g-r color range from -0.96 to 1.65 mag. 
Here, the magnitudes of g and r bands and the color of g-r are all derived from the optical 
surface photometry done by ourselves with elliptical apertures by the SExtractor software 


flBertin &: Arnouts 


19961) . the process of which will be described detailedly in next section. 
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distance 


r-band magnitude 




HI mass 



log[M(HI)] [Mgl 
HI line width 



Fig. 1.— Property distributions of the parent sample of 12423 q;. 40-SDSS DR7 sources. The 
properties are (a)distance with a bin of 10 Mpc, (b)heliocentric radial velocity with a bin 
of 500 km s“^, (c)r-band magnitude with a bin of 0.25 mag, (d)g-r color with a bin of 0.1 
mag, (e)Hl mass with a bin of 0.1 dex, and (f)Hl line width with a bin of 10 km s“^. Here, 
magnitudes and colors in (c) and (d) are derived from our own elliptical-aperture photometry 
described detailedly in Section 3 instead of SDSS photometry . 
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3. Analysis of the optical images of the parent sample 
3.1. sky subtraction 


SDSS analyzes the raw image data, including bias subtraction, sky and PSF 
determination, flat-fielding, and finding and measuring the properties of objects by its 
Photometric Pipeline (PHOTO). For subtraction of the sky background of images, PHOTO 
performs a strategy of a little simplicity. It at first estimates an initial global sky which 
is taken from the median value of every pixel in the image, clipped at 2.32 ct. Secondly, 
PHOTO proceeds to find all the bright objects which are typically those pixels with values 
more than 51 times of the initial global sky estimation and then mask the bright object 
pixels. Third, once the bright objects have been masked, PHOTO determines the same 
clipped median sky value locally within each 256x256 pixel box on a grid with 128 pixels 
spacing and then bilinearly interpolates this sky value to each pixel to construct a sky 
background image. Finally, the sky image is subtracted from the image. 


This strategy of sky background estimation used in the SDSS PHOTO has its 
disadvantage, which has a high risk of considering the large extended outskirts of the 
bright objects as part of the sky background, and thus overestimates the sky background 


of the bright obiec ts, inevitably causing underestimat i on of the luminos i ty and rac 


the bright objects (iLaner et ah 


20671 : 


Lin et ah 


2008 


Hyde et ah 


2009 


He et ah 


ins o 


20131). 


Tis sky subtraction prob 


fjAdelman-McCarthv et al. 


em has also been recognized by the SDSS team themselves 


2006 


20081) . and they quantified the underestimation of the 


brightnesses of galaxies of large angular e xtent due to poor sky subtractio n to exceed 0.2 


mag for galaxies brighter t han r=14 mag flAdelman-McCarthv et al 


proved to be even larger in 


Lauer et al 


2008 ). This bias was 


fj2007|) . which found that for the brightest cluster 


galaxies (BCGs) in their sample of luminous galaxies, the SDSS luminosities and radius are 
strongly biased to low values by excessive sky subtraction. For the BCGs with large effective 
















































14 


radiuses {Re >60 arcsec), the underestimation discrepancies of the SDSS r -band magnhudes 


resulted from the excessive sky subtraction are even larger than 1.0 mag. 


Liu et ah 


fl2008h 


performed accurate sky background subtractions and surface photometry for a carefully 
selected sample of 85 BCGs in the SDSS r band. The comparison of their photometric 
results with those of the SDSS demonstrated that the SDSS pipeline underestimated the 
sizes and luminositie s of BCGs, and th e discrepancies would become larger if the sizes of 


BCGs were larger. 


Hyde et ah 


(120091) compared the SDSS photometric reductions with 


those of their own code for a subset of their full sample of early-type galaxies and expectedly 
found that the SDSS underestimated the sizes and brightness es for large objec ts because of 


the sky subtraction problems of the SDSS PHOTO pipeline. 


He et ah 


( 20131) performed 


accurate sky subtraction and surface photometry for a complete and homogeneous sample 
of bright early-type galaxies (ETCs) and compared their measurements with those from 
the SDSS DR7. They found that the SDSS measurements are on average 0.16 mag and 
maximumly 0.8 mag lower than their own measurements in Petrosian magnitude due to 
the overestimations of the sky background by the SDSS PHOTO pipeline and smaller in 
effective radius. Such underestimations of the luminosities and sizes of the brightest ETGs 
also led to underestimations of the luminosity density and stellar mass density for bright 
ETGs. This recognized issue of the sky background subtraction in the SDSS data releases 
(DRs) from DRl to DR7 has been signihcantly improved by reprocessing all SDSS imaging 
data using a more sophi sticated algorithm for sky background subtracti on in the SDSS 


DR8 ( Aihara et al. 


20111) . However, comparing with the measurements of 


He et al. 


(12Q13I), 


it still underestimates the luminosities of ETGs by about 0.12 mag on average due to the 
overestimations of the sky background subtraction in SDSS DR8. 


Low-luminosity galaxies tend to have low surface brightness than avergae, and 
the flux of such a galaxy with low surface brightness is signihcantly reduced because 
the sky subtraction determination subtracts a substantial fraction of the galaxy light 
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Strauss et al. fj2002h : iBlanton et al. fl2005h . For galaxies with extended low surface 
brightness outskirts which are easily considered as the sky light by the SDSS photometirc 


pipeline, underest imations for 


high as 0.5 mag flLisker et ah 


umino sities are still serious, and the bias can even reach as 


20071) . So, there are limits of the SDSS sky background 


subtraction algorithm, not the SDSS data, but deriving a sample of LSB galaxies from 
photometry alone requires accurate sky background subtraction since an overestimation of 
background like that in the SDSS PHOTO pipeline would bias the number of the true LSB 


galaxies towards a higher value. The phot ometric pipeline of t 


for hnding low surface brightness objects flBlanton et al. 


le SDSS was not optimized 


20051) . Therefore, if we expect 


to derive a reliable LSB galaxy sample from the a.40-SDSS DR7 parent sample by optical 
photometry, we have to carefully estimate the sky background images and do the surface 
photometry ourselves for the SDSS DR7 images of all of the 12423 galaxies in our parent 
sample. 

In preparation, we derive the FITS image hies of the corrected frames (fpC-images) 
with 2048 X 1489 pixels in both g (4686 A) and r (6165 A) bands from the SDSS DR7 
database for all of the 12423 galaxies in our parent sample. The fpC-images are the 
images which have been pre-processed by the SDSS photometric pipeline, including bias 
subtraction, hat-helding, purge of bright stars and corrections for bad pixels (bad columns, 
bleed trails, and those corrupted by cosmic rays ) by interpola t ed va l ues. Then, we ado pted 


a more precise method which has developed by 


Zheng et al. 


fllOOOl) : 


Wu et al. 


(120021) to 


estimate the sky background image for every fpC-image in our parent sample. 


Firstly, we produce a smoothed version of each fpc-image by hltering it with a Gaussian 
function of FWHM = 8 pixels to make the area of each object a bit more extended in the 
image (the smoothed fpC-image). 


Secondly, we use the software SExtractor flBertin &: Arnouts 


19961) to automatically 
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detect all objects with the peak flux higher 2.0a above a global sky background value which 
is simply estimated by the SExtractor itself in the fpC-image. Subsequently, we mask out 
all the objects detected by SExtractor in the fpC-image to produce an object-masking 
image. We then carefully check this object-masking image by eye and unfortunately And 
that the wings of bright stars or the faint stellar haloes of galaxies have not been nicely 
masked. This would absolutely lead to the overestimation for the sky backgroud and then 
underestimation for the luminosity of objects. Through trial and error, we And that if the 
orginal fpC-image is directly used for SExtractor to detect objects, it may be hard to derive 
a nice object-masking image which have the wings of the bright objects or the faint stellar 
haloes of galaxies completely masked. Instead, the smoothed fpC-image generated in last 
step is an optimal choice for SExtractor to produce a nice object-masking image by nicely 
detecting objects. It is worth noting that, empirically through trial and error, the Gaussian 
function with FWHM = 8 seems to be the best choice because it can produce the nicely 
smoothed fpC-image which can be used by SExtractor to nicely detect not only the bright 
inner regions of objects but also the just wings of bright objects and the faint stellar haloes 
of galaxies and then can generate a complete (nice) object-masking file. (e.g. Figure [2] (b) 
is the complete mask file for the orginial fpC-image of Figure [2] (a)). 

Thirdly, we subtracted all of the objects from the fpC-image according to the masked 
areas defined in the complete object-masking file. The object-subtracted image would leave 
us only but sufficient sky pixels from which we could attempt to precisely determine a 
reliable sky background. 

Lastly, we used the object-subtracted image to model the sky background. We 
perform a least-square polynomial fit of low order to the sky pixels of each row of the 
object-subtracted image and replaced the masked pixels of this row by the fitted values. 
Here, the reason why we restrict the fits to low-order polynomials is that it can avoid 
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introducing spurious fluctuations to the masked regions in the fitted sky by interpolations. 
This fitting process was performed firstly row by row and then column by column. The 
individually derived row-fitted and column-fitted sky images were then averaged, and this 
averaged image was later smoothed with a median Altering box 31 x 31 pixels in size to 
eliminate any small artifacts from the modeling process. This smoothed sky images is 
Anally adopted as the sky background (Figure [2] (c)) and is subtracted from the original 
fpC-image (Figure |2] (d) is the sky-subtracted image of the original fpC-image shown as 
Figure [2] (a)). Note that the sky backgroud (Figure [2] (c)) shows a gradient across the 
frame as expected. Comparing with a straight 2D background fit which has been fou nd 


to systematically underflt or overflt certain regions of the images fjZheng et ah 


199911 . 


the method of fitting the sky piecewise in a row-by-row and column-by-column fashion 
using low-order polynomials can ensure us to predict the sky background underneath the 
object-masked regions in a mutually orthogonal manner, which not only permits reasonable 
interpolations of sky und er the galaxy region but also fit the inherent lumpiness of the sky 


at low surface brightness (jZheng et ah 


19991 ). Besides, this bidirectional A 


alreac 


2002 


y been checked a nd success 


Lin et ah 


Liu et ah 


2008 


2003 


Liu et ah 


Chonis et ah 


ully used in many 


2005 


2011 


Duan 


2006 


Mao et ah 


papers fIZheng et a 


Cao fc Wu 


2007 


ting method has 


1999 


Wu et ah 


Li et ah 


2007 


20i4 . 


Following the speciflc steps listed above, we have accurately estimated sky background 
maps and then subtracted the sky backgrounds from the fpC-images for all galaxies in our 
parent sample in both g and r bands. To check the quality of our sky subtraction, we have 
made necessary tests. As an example, for the galaxy presented in Figure |21 Figure |3] shows 
distributions of counts in the sky-subtracted frame (the solid black) for all unmasked pixels 
of the whole frame (global; the left-hand panel) and of the local vicinity around our galaxy 
(local; the right-hand panel). The around local vicinity is deflned by the region between 
boundaries of two square boxes repectively sized 250 x 250 pixels and 500 x 500 pixels from 
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the galaxy center. If our sky background model is successful, the counts are expected to 
follow a Gaussian distribution with a mean close to zero both globally and locally. This is 
indeed the case, as can be seen from the distribution represented in solid black in Figure El 
The count distributions for the whole frame and the local region in the SDSS r band are 
both well htted by the Gaussian functions with the mean values very close to 0 ADU. The 
dispersions of the Gaussian distributions are respectively 6.83 ADU for the whole frame 
and 6.74 ADU for the local region. Additionally, the distributions represented in dashed 
grey in Figure [3] are count distributions in the frame without sky background subtraction 
by our accurate row-by-row and column-by-column htting method but only with the simple 
global or local mean value subtraction for all unmasked pixels of the whole frame and of the 
local region. These dashed grey distributions are plotted for a further demonstration of the 
nicety of our sky subtraction. It is clear that the mean values of the Gaussion distributions 
for the frame with our sky background subtraction (the solid black) are much closer to zero 
than for the frame with no accurate sky background subtraction but only with the simple 
global mean value subtraction (the dashed grey) both globally and locally (0 vs. -1.63 ADU 
for the whole frame and 0.03 vs. 0.09 ADU for the local region). Evidently, this row-by-row 
and column-by-column method of sky background subtraction can well recover the sky 
background model, but the most noticable problem of this method, comparing to a straight 
2D htting, might be a little more time-consuming as a result of its relatively complicated 
htting algorithm. 


3.2. Photometry by SExtractor 


I n the system o 


used fIStrauss et ah 


the SDSS Petrosian photometry, the Petrosian circular apertures are 


2 OO 2 I) . Although the circular aperture can work well for objects with 


small angular extent or spherical-like morphology, it is not an optimal choice for galaxies 
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(c) sky background 


Fig. 2.— Sky subtraction process for a galaxy (AGC 4130) in our parent sample as an 
example. The objective galaxy AGC 4130, as the edge-on galaxy in this frame locates 
at the bottom-right region of every picture. The four pictures respectively show (a)the 
original SDSS DR7 fpC-image in r band for the galaxy, (b)the object-masking hie, (c) the 
sky background map from htting, and (d)the hnal sky-subtracted image . 
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Fig. 3.— Distributions of count for all of the unmasked pixels in the whole frame (global; 
the left panel) and in the local vicinity region around the galaxy (local; the right panel). 
This local vicinity region is dehned as the image area between boundaries of two square 
boxes sized respectively 250x250 pixels and 501x501 pixels from the galaxy center. The 
solid black represents distributions for the skysubtracted frame by our piecewise row-by-row 
and column-by-column htting method, and, for a further comparison, the dashed grey are 
for the frame without our accurate sky background subtraction but only with the simple 
global or local mean value subtraction for all unmasked pixels of the whole frame and of the 
local region. . 




21 


such as most of the galaxies in our parent sample shown in Figure IH which usually have 
large angular extent, irregular morphology or edge-on shape. For objects like those shown 
in Figure 01 the circular apertures are too small in size to include all the inherent intrincic 
light from the objects that we are interested in, or for objects shown in Figure IU the circular 
apertures are so large that they inevitably involve the light from the adjacent objects. 

As large numbers of LSB galaxies are morphologically similar to late type (Sc and later) 


spiral galaxies with amorpho us or fragmentary and 


similar to irregular galaxies flMcGaugh et al 


hint spiral patterns, or morphologically 


I 995 I) . Therefore, the circular apertures 


are obviously not the best choices for LSB galaxies. In this paper, we will do surface 


photo metry with the elliptical apertures by the SExtractor software flBertin fc Arnouts 


19961 1 for the sky-subtracted frames of every galaxy in our parent sample. 


The SExtrator package is a source-detection and photometry routine. For photometry, 
there are six different approaches (isophotal, corrected-isophotal, automatic, best, hxed- 
aperture, and Petrosian) in the SExtractor software. Thereinto, the automatic ap erture 


magni tudes (AUTO), inspired by Kron’s “ hrst moment” algorithm (see details in 


Kron 


(1l980l) j. are intended to give the most precise estimate of “total magnitudes”, at least for 
galaxies. SExtractor’s automatic aperture is an flexible and accurate elliptical aperture 
whose elongation e and position angle 6 are dehned by the second order moments of 
the object’s light distribution. Then, within this aperture, the characteristic_radiusri is 


dehned as tha 


Infante 


weighted by the light distribution function (ri= ^^^^) 

(119871 ) have verihed that for stars and galaxy prohles convolved with Gaussian 


Kron 


fll98n[l and 


seeing, more than 90% of the hux is expected to lie within a circular aperture of radius 
kri if k=2, almost independently of their magnitudes. This picture remains unchanged 
if they consider an ellipse with ekri and ^ as the principal axes. By choosing a larger 
k=2.5, more than 96% of the hux is captured within the elliptical aperture. So, we 
keep k=2.5 that is also the default setting of SExtractor during the automatic elliptical 
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aperture photometry flBertin fc Arnouts Ill996l ). More details about the Kron ra. d ius an d 


the automatic aperture photometry used in SExtractor can be derived in 


Infante 


dlMZl); 


Bertin &: Arnouts 


Kron 


fll996l 


(1198(11 : 


Therefore, we would perform the photometry by SExtractor for the sky-subtracted 
fpC-images in both g and r bands of every galaxy in our parent sample. Firstly, the 
sky-subtracted fpC-images i n r band of all the galax ies in our parent sample were fed 


into the SExtractor routine (jBertin &: Arnouts 


19961) to derive the r-band magnitude and 


galaxy centeral coordinate of each galaxy. Then, we used the same galaxy center and 
aperture as those of the r-band image to measure the g-band sky-subtracted fpC-images 
and derived the g-band magnitude of each galaxy. We show examples of the SExtractor 
AUTO elliptical aperture at the SExtractor-determined galaxy center (white ellipses in 
Figure H, comparing with the SDSS circular Petrosian aperture at the SDSS-determined 
galaxy center (red circles in Figure H]) for 4 galaxies in our sample. Obviously, the AUTO 
photometry appears indeed to be more appropriate for those galaxies because the AUTO 
ellipses could capture all of the light from the objective sources themselves but exclude the 
light from other nearby objects. However, the SDSS Petrosian circles adopted in SDSS 
Petrosian magnitudes seem to have excluded amounts of light from the objective sources 
themselves or even included light from ajacent sources if the objective galaxy was to be 
wholely included in a circular aperture. 


The photometry for all the galaxies using SExtractor was automatically done in a 
batch mode, and then we made a visual check on the AUTO aperture photometric results 
and redid the AUTO photometry by correcting the inappropriate apertures. Finally, the 
resultant g band magnitude and g-r color are previously shown in Figure [T] (c) and (d). 
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Fig. 4.— Comparisons between AUTO elliptical apertures used by SExtractor and the 
Petrosian circular aperture used by SDSS in the r band for galaxies from our parent sample. 
The white ellipses are the SExtractor elliptical AUTO apertures which are centering at 
the SExtractor-determined galaxy center. For comparisons, the SDSS cicular Petrosian 
apertures (red circles) which are centering at the SDSS-determined galaxy centers have been 
also plotted . 
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Fig. 5.— Examples of SDSS petrosian aperture (SDSS Petro) and the Kron elliptical aper¬ 
tures (SEx AUTO) from SExtractor for 2 galaxies from our parent sample. The SDSS Petro 
circular aperture are shown as the red circles, and the SEx AUTO elliptical apertures by 
the SExtractor routine are shown as the middle white ellipses. Obviously, the SEx AUTO 
elliptical apertures seem to be much more approriate for the LSB galaxies . 
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3.3. Geometry by Galfit 


Although, besides photometry, some interesting fitted parameters can also be made 
by SExtractor currently, these fitted results are only the rough guess and estimate of the 
parameters because the advantage of SExtractor is its photometry instead of geometric 
fitting by which SExtractor only aims to give a quic k look for the int eresting fitted 


parameters. So, we choose to use the Galfit software flPeng et al. 


2 OO 2 I) . which is good at 


galaxy fitting. By setting a single radial profile function or a combination of a number of 
functions, e.g. the Sersic, exponential, Nuker models and others which are allowed in most 
literatures and artifically setting the initial values for a set of input parameters, Galfit starts 
a nonlinear least squares fitting. During the fitting, minimization of residuals between 
model and image is done by using the Levenberg-Marquardt downhill-gradient method. 
The process of minimization iterates until convergence is achieved. Then, the solution in 
the case of convergence should be regarded as the optimum solution in the parameter space. 


Therefore, after surface photometry by the SExtractor software, all of the sky- 


subtracted images in g and r bands of t 

le 12423 ga 

axies of our parent sample were fed 

into the galaxy-fitting procedure Ga 

of LSB galaxies lack strong bulges ( 

[fit ( 

Peng et al. 

2002 

) in a batch mode. A 

lS maiorities 

de Blok et al. 

1995; 

Beiiersbergen et al. 

1999 

), a 


decomposition into bulge and disc is not essential for LSB galaxies, especially for the 


more disk-dominated LSB galaxies which are preferred in our future science goal. So, we 
fit each galaxy in our parent sample only with an exponential profile function, and set 
the parametric results from the SExtractor to be the initial values for the set of input 
parameters (galaxy magnitude, disk scale length, b/a ratio and position angle) for Galfit, 
which would finally determine the best-fit values for the interesting parameters including 
the disk scale length in pixel(a), axis ratio (g), and inclination angle {i) for each galaxy of 
our parent sample and generate a triplet of thumbnails including the original galaxy, the 
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exponential model fit and the residual image (see Figure |6]). 


Although the best-£t results are derived in the case of convergence and minimization, 
problems with local minima and numeri cal degeneracies b etween some parameters may be 


present as the parameter space is large (jPeng et ah 


20021). To reinforce our conhdence in 


the parameter optima, we use Monte Carlo simulations to search the parameter space for 
each galaxy of our parent sample. We make 200 Monte Carlo simulations for each galaxy of 
our parent sample by the following steps. 

1. Set the exponential function as the galaxy radial prohle function for Galht. 

2. Randomly select a set of parameters (magnitude, disk scale length, b/a ratio, position 
angle), drawing from distributions of these parameter values centered on the best-£t values 
previously determined by Galht htting. 

3. Using these randomly-selected parameters as initial guesses for the set of input 
parameters, minimize the x^ until the convergence is reached by Galht. 

4. Repeat steps 1-3 for 200 times to see whether they return the same optimized solution as 
that derived by using the SExtractor results to be the initial guesses for Galht or to other 
equally plausible ones. 


For convenience, we use the term of “best-ht solution” to denote the best-ht solution 
from Galht by using the SExtractor results as the initial guesses, and the term of “simulated 
solutions” to denote the solution from Monte-Garlo simulations by Galht which use 
randomly-selected parameters as the initial guesses. We compare the “best-ht solution” 
with the “simulated solution” for galaxies of our parent sample and hnd good agreement 
between the two solutions for every galaxy, which strongly verihes the validity of our best-ht 
solutions for the useful galaxy parameters. For an example, in Figure [71 we take the same 
galaxy (AGG 572) as that shown in Figure [6]which shows the galaxy-htting process of Galht 
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to present the comparison between the “best-£t” and the “simnlated” solntions. In Fignre 
[3 comparisons between the “best-£t solntion” (the solid line) and the “simnlated solntions” 
from 200 Monte-Carlo simnlations (the dots) are shown for the two usefnl parameters in 
this paper of the disk scale length and axis ratio respectively from top to bottom. For each 
panel from this figure, the “simulated solutions” from 200 Monte-Carlo simulations strongly 
converge at their mean value (denoted as “Simulation mean and sigma” in Figure [7]) and 
the mean value of the “simulated solutions” is well consistent with the “best-£t” solution 
from Galht (denoted as “Galfit” in Figure [7|)- It is worth noting that in the bottom 
panel of Figure 7, it looks like that the axis ratio b/a values from simulations 
are 0.8108, 0.8109, and 0.8110 only, with nothing in between. Acturally, it is 
because that the Galfit software has rounded its final output results to only 
four decimal places, so any difference between the b/a values after four decimal 
places can not be apparent to us. This further strengthens the consistency of 
the b/a values from 200 Monte-Carlo simulations as the differences are at least 
after three decimal places. 


3.4. Central Surface Brightness 


The typically used photometric parameter to separate the 

brightness regime of galaxies is the central surface brightness o 

high and lovv 

' the disc in t 

surface 

le B band. 

unlBl (Freeman 1970: Imnev et ah 1996: 

O’Neil et ah 1997a: 

Zhong et ah 

2008). So in 

this subsection, we will calculate /xo(B) for each galaxy in our parent sample. Deduced from 


the exponential function of the radial profile of the disk component of a galaxy, the central 
disk surface brightness can be expressed as eq.(la) where po refers to the central disk surface 
brightness in mag arcsec“^, m refers to the apparent magnitude in mag, a represents the 
disk scale length in pixel and A is the area of one pixel. For the SDSS image, A equals to 
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0.396x0.396 in arcs ec^/pixel. Following 


O’Neil et al. 


Zhong et al. 


fll997ah : 


Trachternach et al. 


fl2006h : 


(120081) ■ the central snrface brightness shonld fnrthermore be corrected by 


inclination and cosmological dimming effects, so it can be hnally expressed as eq.(lb), 
where q and z refer to the axis ratio and redshift, respectively. 


/io = m + 2.5 \ogiQ{2TTa^ A) (la) 

/To = + 2.5 logio(27ra^v4g) - 10 logio(l + z) (lb) 

ho(5) = /io(5') + 0.47(/io(5') - + 0.17 (Ic) 


Using the magnitnde measnred with the AUTO apertnre by SExtractor and 
snbseqnently Galactic-extinction corrected by ns as m, the disk scale length in pixel (a), 
the axis ratio (g) from the exponential profile fittings by Galfit, the redshift (z) dednced 
from the Hi velocity in the a.40 parametric catalog, and the pixel area of A= 0.396x0.396 
in arcsec^/pixel for the SDSS image, we firstly calcnlated the central disk snrface brightness 
in both g and r bands, /io(g) and poij), in nnits of mag arcsec“^ according to Eq.(lb). 
Then, the B-band central disk snrface brightness was calcnlated from po(s) and po(^) by 
nsing (eq.lc), which st ems from the f ilter t ransformation relation between g, r and B band 


(Eq.(2a)) derived from 


Smith et al. 


( 120021 ). 


4. Low Surface Brightness Galaxies 
4.1. LSB galaxy sample 

LSB galaxies are commonly defined as a popnlation of diffnse galaxies whose central 
snrface brightness in B band, /io(B), falls below a specific threshold valne, which is 
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21.65±0.30 in 


Freeman 

(1970 

), 22.0 in 

O’Neil et ah 

n997a 

), and 23.0 in 


Impev fc Bothun 


fjl997n . In general, the most common threshold values found in the literature are between 
21.5 and 23.0 mag arcsec^. In this paper, we adopted a threshold of > 22.5 mag 

arcsec“^ for the B-band central surface brightness and put additional constraints on 
the axis ratio (b/a > 0.3) to select a sample of LSB galaxies with no edge-on galaxies. 
With these two constraints, we hnally constructed a sample of LSB galaxies consisting 
of 1129 non-edge-on LSB galaxies from the a.40-SDSS DR7 survey. This sample is a 
relatively unbiased sample of the disc-dominated LSB galaxies and it is complete both in 
Hi observation and in optical magnitude within the magnitude limit of SDSS photometric 
observations. 


4.2. Statistical properties of LSB galaxy sample 


4.2.1. Optical properties 


Distributions of the B-band central surface brightness (/io(B)), absolute magnitude 
(M(B)), the B-V color (B-V), and the scatter of /io(B) against B-band apparent magnitude 
(B) for all galaxies in our LSB galaxy sample are shown in Figure [8] (a)~(d). 

The B-band central surface brightness of all these LSB galaxies (Figure [H] (a)) are 
distributing from 22.5 to 28.3 magarcsec~‘^, with 53% (597/1129) betweeen 22.5 and 23.0 
magarcsec~‘^, 34% (388/1129) between 23.0 and 24.0 magarcsec~‘^, 9% (103/1129) between 
24.0 and 25.0 magarcsec~‘^, and the left 4% (41/1129) with the B-band central surface 
brightness fainter than 25.0 magarcsec~‘^. 


The coverage of the B-band absolute magnitude (Figure [8] (b)) of all the galaxies in 
the LSB galaxy sample are from -27.0 to -12.3 mag. There are 43 galaxies in ou r LSB 


galaxy sample being instrinsically very faint (M(B)>-17.3 mag] 


Poggianti et ah (1200llH 
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We present the distribution of the B-V color of the LSB galaxy sample in Figure lUj^c), 
which reveals that our sample covers from -0.22 to 1.88 mag in B-V color and it has 98.4% 
galaxies bluer than B-V=0.75 mag. This strongly probes that this sample is a bluish 
sample of LSB galaxies. Additionally, we depicted the distribution of B-band central surface 
brightness againt apparent magnitude in Figure IH] (d), from which we can see a trend that 
fainter galaxies tend to have fainter central surface brightness. 


4-2.2. 21cm H\ properties 


We derived the radio information of all the galaxies in our LSB galaxy sample from the 
ALFALFA catalogue and showed the distributions of them in Figure IH] (e)~(h). 

As shown in Figure | 8 ] (e), the common logarithm of Hl-gas mass (log M{Hi)/M q) of 
the LSB galaxy sample ranges from 6.11 to 10.36 dex, with 32% (365/1129) having high 
mass of Hi gas (log M(Hl)/M© > 9.5 dex), 63% (712/1129) having medium mass of Hi 


gas (log M(Hl)/M 0 - 7.7~9.5 dex), a nd the left 5% 


of Hi gas (log M(Hl)/M© < 7.7 dex; 


Huang et ah 


(52/ 


129) galaxies having low mass 


(1201211 1 . Such a distribition indicates 


that LSB galaxies are mostly rich in gas with high or medium mass of Hi. The velocity 


width of Hi line (W50) covers from 11.0 to 443.0 kms ^ (Figure M (f)), 


40% 


450/1129) have narrow Hi line velocity width (W 50 < 80 kms 


-1. 


rom which on ly 


Huang et ah 


(120121) 1. The distribution of ratios of Hi mass and B-band luminosity (measured by our own 
photometry) are depicted in Figure | 8 ](g). Our LSB galaxy sample has a median M(Hl)/Ls 
= 0.87, which suggests that most LSB galaxies of this sample are gas-rich. In Figure [ 8 ] (h). 


we show the distributions of the B-band central surface brightness against Hi mass. 
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4-2.3. Environmental properties 


Evironment is a critical factor which affects the star formation and evolution of 
galaxies. In this subsection, we made a simple study on the environment where galaxies in 
our LSB galaxy sample live. 

In order to have a overview, we provide graphic illustrations of distribu tions of our 


LSB g alaxies overplotted on the large-scale structure plot (Figures 5 and 6 in 


Havnes et ah 


(120111 )) of the ALFALFA survey in the local universe. As an example. Figure [9] in our 
paper shows cone diagrams of a four-degree wide slice of the LSB and a.40 galaxies in 
the ALFALFA spring and fall sky centered on decl.= -|-26°, including the full ALFALFA 
bandpass redshift range cz < 18000 km s“^. Blue open circles mark the locations of galaxies 
detected by ALFALFA survey, while red open circles denote galaxies in our LSB galaxy 
sample. 


To study the environment of galaxies i n our LSB galaxy 


we follow the method of Ending clusters in 


Wen et ah 


sample quantitatively. 


(120091) . For each galaxy in 


our LSB galaxy sample at a given redshift, z, we counted the number of all galaxies 
detected by the photometric observation of SDSS DR7 within a radius of 1.0 Mpc from 
the center represented by the galaxy postion and a photometric redshift gap between 
z ± 0.04(1 -|- z). Fortunately, the SDSS provides two useful functions for users. One is 
the f Get Nearby Obj ARE q function which returns a table of all objects detected by SDSS 
within a radius in arcmin of a given equatorial point, and the other is the fCosmoDa 
function which returns the angular diameter distance at a given redshift. So with the help 
with these functions, we can easily derive the number of all objects within that radius and 
redshift space detected by SDSS photometric survey. We show the simple environmental 
properties of our LSB galaxy sample in Figure [TOl which shows the distributions of the 
number counts of galaxies within a radius of 0.5 Mpc in panel (a), 1.0 Mpc in panel 
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(b), and 3.0 Mpc in panel (c) and scatter distribntion of the detectable galaxy nnmber 
connt with a radins of 1.0 Mpc against the Hi mass for this LSB galaxy sample in panel 
(d).Obviously in Figure ITUrbh about 92% out of the total galaxies in our LSB galaxy 
sample have less than 8 galaxies around within the 1.0 Mpc radius and z ± 0.04(1 + z) 
redshift space. It strongly distinguishes from the co unt distribut i on of the member galaxy 


candidates of clusters within 1.0 Mpc in Figure 4 in 


Wen et ah 


(120091) . which distributes 


totally from greater than 8 (grey dashed line) to 50 with the peak at 16 (grey solid line), 
comparing to the peak of 2~3 for our LSB galaxy sample. Even if the radius is increased 
to 3.0 Mpc, the portion of galaxies having less than 8 detectable neighbours is still as high 
as 65%. This dehnitely indicates that LSB galaxies are more likely to reside in the held 
environment. Additionally, from Figure [T0f a).fb) and (c), the parent sample has nearly 
the same number of neighbours as the LSB sample. This is large because that the parent 
sample itself is a gas-rich galaxy sample from the ALFALFA Hi survey as galaxies rich in 
gas are generally favors low-density environment. In Figure [TOl i. we show distribution of 
the neighbouring galaxy counts againt the Hi mass for our LSB galaxy sample as black 
dots. As mentioned before, the Hi mass of our LSB galaxy sample ranges from 6.5 to 11.0 
dex in log M(Hl)/M©. For a clear understanding of the probable trend between the count 
of neighbouring galaxies and the Hi mass, we divided the Hi mass range of the LSB galaxy 
sample into 9 bins with a binsize of 0.5 dex in log M(Hl)/M©, and investigated the relation 
(the grey broken line in Figure [TOtH between the mean Hi masses and mean number counts 
of neighbouring galaxies of the 9 mass bins, which does not give a visible relation between 
the mass and mean counts of neighbouring galaxies, but shows that the mean counts of 
neighbouring galaxies in all massbins appear lower than 5. 
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5. Discussion: impact of bulges on fitting results 


As LSB galaxies are assumed to be disc-dominated and lack strong bulges, we only 
adopt a single exponential function as the model in Section 3 to perform a fast search of 
LSB galaxies from our parent sample by Galfit. This fast method of using only a single 


expon ential profile during the Galfit fitting has already been used by 


Trachternach et ah 


(120061 ) to search LSB galaxies from gas-rich galaxies of a blind optical survey. This fast 


method tends to be good at detection of the disc-dominated LSB galaxies. However, it 
may preclude detection of the bulge-dominated LSB galaxies. In Figure [TT], we show the 
distribution of the bulge-to-total ratio (B/T) for 1129 galaxies of our LSB galaxy sample 
(dashed black), and for a comparison the B/T distribution of all of the 12423 galaxies of 
our parent sample is overplotted (solid grey). Both histograms are normalized to have 
a maximum value equal to 1.0. It is worth noting that these B/T values are derived 
from the parameter ’’fracDev^” from the SDSS DR7 catalogue, which is defined as the 
ratio of luminosity contributed by the bulge relative to the total luminosity of galaxies 
in r band. As is once proposed that B/T = 0.4 can be roughly considered as a diyision 


between late-ty pe (disc) and early-type (E/SO) galaxies (jSimien fc de Vauconlenrs 


1986 


Li et ah 


20071) . from Figure HH it can be concluded that among the 1129 galaxies of 


our LSB galaxy sample, 50.6% (571 galaxies) are pure disc galaxies (B/T =0), 43.( 

(492 galaxies) are disc-dominated galaxies (0<B/T<0.4) and only 5.8% (66 galaxies) are 
bulge-dominated galaxies (0.4<B/T<1.0). Gomparing with the composition ratio of our 
parent sample, 26.3%, 50.2% and 23.5% of which are respectively pure disc, disc-dominated 
and bulge-dominated galaxies, our LSB galaxy sample is more of a disc-dominated LSB 
galaxy sample that is deficient in the bulge-dominated LSB galaxies. 


In this disc-dominated LSB galaxy sample, it is correct to fit the radial profile of the 
pure disc galaxies with only a single exponential profile. However, for the galaxies that do 
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have bulges in this LSB galaxy sample, the derived central disk surface brightnesses might 
be biased towards brighter values by htting only with a single exponential profile. Instead, 
a two component (Sersic + exponential) study would be much better for these LSB galaxies 
with bulges. To carry out a test, we randomly take 20 galaxies with different levels of 
bulges (0.05 < B/T <1.0) and morphologically apparent disks from our LSB galaxy sample 
as the test sample. Then, for this test sample, we fit the galaxies one by one by Galfit with 
a combine of two components which are the disk component expressed by an exponential 
radial prohle and the bulge component expressed by a Sersic radial prohle. Comparing with 
the previous results of htting these galaxies by Galht with only a single exponential radial 
prohle, the two component htting results seem to be better, especially for the galaxies 
with morphologically larger bulges because the bulges of these galaxies can be subtracted 
more cleanly as can be seen from the residual images of the two component htting results. 
In Figure [121 we show images of both the single component htting (top panel) and two 
component htting (bottom panel) by the Galht for one galaxy with B/T equal to 0.61 
from the test sample for an example. As can be obviously seen from Figure [121 fh® fwo 
component study can ht the bulge of the galaxy more sufficiently than the htting only by a 
single exponential prohle. 

Being ahected by the existing central bulge, the central disk surface brightness would 
be biased to a brighter value if only the single exponential prohle is used to ht this galaxy. 
The top panel of Figure [131 shows comparisons between the derived central disk surface 
brightnesses for the 20 galaxies of the test sample respectively from the single exponential 
prohle httting and the two component study by Galht, which shows that the central disk 
surface brightnesses are generally overestimated (to brighter values) by htting only with a 
single exponential prohle. The bottom panel of Figure [131 shows that the overestimations 
of central disk surface brightnesses are becoming higher for galaxies with relatively larger 
bulges. 
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Our LSB galaxy sample is established basing on the central disk surface brightness 
derived from the Galht htting by only a single exponential prohle, so the central disk 
surface brightnesses of this LSB galaxy sample are overestimated (to brighter values) 
more or less, especially for those galaxies with large central bulges or light concentrations. 

The stronger the bulges, and the larger the overestimations are (Figure I13jl . 
Nonetheless, galaxies with stronger bulges or light concentrations are not 
dominant in our LSB galaxy sample (Figure 11111 . Quantitatively, from the 
bottom panel of Figure I13L for LSB galaxies with weak or no bulges (B/T < 
0.2), which contributes 84.7% (956/1129) to the sample size, the impact of 
introducing a bulge component on the derived central surface brightness is less 
than 0.2 mag arcsec”^, although the effect can be larger than 0.2 mag arcsec”^ 
for the left 15.3% fraction (173/1129) of galaxies in our sample when bulges are 
stronger (0.2<B/T < 1). 

In the case of overestimation, the true central disk surface brightnesses for those galaxies 
already in our LSB galaxy sample should be even fainer than currently, so these galaxies 
currently in our sample should still be members of our LSB galaxy sample. Although 
LSB galaxies are rare to have large bulges or central light concentrations, 
such bulge-dominated LSB galaxies do exist. However, for such galaxies 
which should essentially be true LSB galaxies with relatively large bulges or 
central light concentrations in the parent sample, they tend to be mistakenly 
identified as non-LSB galaxies only due to the overestimation of the central 
disk surface brightness (to brighter values than the LSB galaxy threshold ) by 
our single exponential profile fitting method instead of a two-component fitting. 
Therefore, we have to acknowledge that our LSB galaxy sample in this paper is more of 
a disk-dominated LSB galaxy sample and it might lack the true LSB galaxies which have 
strong bulges or large central light concentrations as a result of the overestimation of central 
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disk surface brightness by using only a single exponential component fitting. This can be 
evidenced by the large contribution (84.7% ) of disk-dominated galaxies (B/T 
< 0.2) to our LSB galaxy sample. 

Although the htting method of using only a single exponential prohle tends to 
preclude the detection of LSB galaxies which have relatively large bulges or central light 
concentrations, it is fast and valid for the fast majority of the LSB galaxies because LSB 
galaxies are believed to be mostly disk-dominated galaxies. In our current LSB 
galaxy sample, the disk-dominated galaxies (B/T < 0.2) contribute up to 84.7% 
to the sample size. Additionally, seen from Figure IllL the B/Ts of this LSB 
galaxy sample distribute essentially the same way as those of all the parent 
sample do at the small B/T end. This consistency may also strengthen the 
validity of this fast method to select the disk-dominated LSB galaxies from all 
of the disk-dominated galaxies in the parent sample as completely as possible. 

In view of the above-mentioned fact that our current LSB galaxy sample is a 
disk-dominated sample and might be dehcient in the true bulge-dominated LSB galaxies, 
in the future, we are planning to make a two-component study for galaxies with bulges 
or central light concentrations in our parent sample and calculate their central disk 
surface brightnesses with as possible as small bias. Basing on the new central disk surface 
brightnesses which are expected to be systematically fainter than the ones derived by only 
a single exponential prohle study, we believe that more true LSB galaxies with bulges or 
central light concentrations would be picked out and join the LSB galaxy sample. Then, we 
will sophisticatedly compare our current LSB galaxy sample with the future sample derived 
by the two-component study, for which it is believed that the contents of Section 4.1 and 4.2 
in this paper will change. However, these two sections now could well present the features 
of our current LSB galaxy sample which is positioned as a non-edge-on disk-dominated LSB 
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galaxy sample selected from the a.40-SDSS DR7 survey. As ALFALFA is a Hi survey, our 
parent sample should be itself composed of more gas-rich galaxies, so this would surely lead 
our current LSB galaxy sample towards a sample of having more gas-rich LSB galaxies and 
lacking sufficient gas-poor LSB galaxies. 


6. Summary 


LSB galaxies are objects with central surface brightnesses at least one magnitude 
fainter than the night sky. They are galaxies of a ty pe that encompasses many of the 


■‘extremes” in galaxy properties (jO’Neil et ah 


20041 1 so they play an important role in 


understanding the galaxy formation and evolution of the uninverse. As most LSB galaxies 
are rich in gas, the surveys of gas-rich galaxies would be good samples for us to efficiently 
select LSB galaxies. 


The ALFALFA survey is a blind extragalactic 21cm Hi survey in the local universe (z 
< 0.06) for 7000 deg^ sky area, 40% of which has already been released to the public as the 
a.40 catelogue. So the a.40-SDSS DR7 sample consisting of 12423 SDSS-primary galaxies 
provides us one of the best survey combinations to select a relatively unbiased sample of 
LSB galaxies which are mostly gas-rich in the local universe. Therefore, we expected to 
dehne a LSB galaxy sample from the a.dO-SDSS DR7 sample in this paper. 


The selection of LSB galaxies is sensitive to the sky backgrounds. As sky backgrounds 
are systematically overestimated for galaxies by the SDSS photometric pipeline, especially 
for those bright galaxies or galaxies with extended low surface brightness outskirts, we 
accurately re-estimated the sky background of SDSS images in both g and r bands 
for each galaxy of the a.40-SDSS sample, using a careful method of row-by-row and 
column-by-column htting. After sky subtraction, we did the surface photometry by the 
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SExtractor software and fitted the geometry by the Galfit software for each galaxy in g and 
r bands. Basing on both the photometric and geometric resnlts derived by onrselves, the 
central snrface brightness in B band (/io(B)) conld be calcnlated for every galaxy. Then we 
selected galaxies with /io(B) > 22.5 mag arcsec“^ and the axis ratio b/a > 0.3 to establish 
a sample of LSB galaxies from the a.dO-SDSS DR7 snrvey. 


This LSB galaxy sample, consisting 1129 galaxies from the a.40-SDSS DR7 sample, is 
a relatively nnbiased sample of LSB galaxies and it is complete in both the Hi observation 
and the optical magnitnde within the limit of SDSS photometric snrvey. This LSB galaxy 
sample spans from 22.5 to 30.0 for /io(B) with a large fraction (87%) in 22.5 < /io(B) < 24.0 


mag arcsec”^, 9% in 24.0 < /xo(B) < 25.0 mag arcsec“^, and the left 4% fainter than 25.0 
mag arcsec”^, from -27.0 to -12.5 mag for the absolnte magnitnde in B band (M(B)) with 
43 faintest galaxies (M(B) > -17.3 mag), from -0.2 to 1.9 mag for the B-V color with 98.4% 
blner than B-V=0.75 mag indicating that it is a bine LSB sample. In the aspect of Hi 
properties, a large portion (63%) of the sample has high mass of Hi (M(Hl) > 10^'^M©) and 
only a few (1%) has very low mass of Hi (M(Hl) < 10^ '^M©). These indicate that members 
of this LSB galaxy sample are mostly gas-rich, with a median M/^j/Ls >1. Additionally, we 
investigated the environment of galaxies in this LSB galaxy sample. We connted the nnmber 
of all neighbonrs of the central galaxy detected by the SDSS photometric snrvey within a 
radins of 0.5Mpc, 1.0 Mpc and 3Mpc and a photometric redshift gap between z ± 0.04(1 
-|- z) and made a statistic on the nnmber connts of neighbonrs. The distribntions of the 
neighbonr connts strongly evidenced that LSB galaxies prefer to reside in the environment 
of low density, comp aring with the neighbonring galaxy connts of clnster candidates shown 


m 


Wenetah 


(120091 1. 


However, this LSB galaxy sample has its own drawback. As discnssed in Section 5, 
the selection of this LSB galaxy sample is basing on the central disk snrface brightnesses. 
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which are derived by performing only a single exponential prohle £t to the galaxies of the 
q;. 40-SDSS sample. This fitting method of only using a single exponential model is fast 
and valid for the fast majority of the LSB galaxies. However, it indeed overestimates the 
central disk surface brightnesses for galaxies with bulges or central concentrations of light 
by comparing with the results of a two-component study for these galaxies in our LSB 
galaxy sample. For LSB galaxies with weak or no bulges (B/T < 0.2), which 
contributes 84.7% (956/1129) to the sample size, the impact of introducing 
a second component (Sersic) on the derived central surface brightness is less 
than 0.2 mag arcsec“2, although the effect can be larger than for the left 
15.3% fraction (173/1129) of galaxies in our sample which have stronger bulges 
(0.2<B/T < 1). The stronger the bulges, and the larger the effect will be. 
Nonetheless, fraction of galaxies with very strong bulges (5.8% with B/T=l)is 
very small in our LSB galaxy sample. As the present central disk surface brightnesses 
systematically bias towards brightner values and should be fainter than currently after a 
two-component study for galaxies with bulges in our LSB galaxy sample, our LSB galaxy 
sample is a sample consisting of more disk-dominated galaxies but being deficient in the 
true bulge-dominated LSB galaxies which are rare but do exist. In a word, the LSB galaxy 
sample defined in this paper is a relatively unbiased sample of gas-rich, non-edge, and 
disk-dominated LSB galaxies from the overlap between the ALFALFA 21cm Hi survey and 
the SDSS DR7. 

We would like to thank the referee for the helpful suggestions. We also thank ALFALFA 
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Fig. 6.— Exponential model fit to the galaxy AGC 572 of our parent sample by the Galht 
routine. The left, middle, and right panels are respectively the original galaxy image (SDSS 
g-band), the exponential prohle model, and the residual image.. 
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Fig. 7.— Comparison between the “best-fit solution” and the “simulated solutions” from the 
Galht routine for the galaxy AGC 572. From top to bottom is respectively the comparison 
for the useful parameters of disk scale length (rs) and axis ratio (b/a). The “best-£t solution” 
is represented by the solid line in each panel and the dots are representing the “simulated 
solutions” from 200 Monte-Carlo simulations. In the bottom panel, b/a values seem to be 
only three discrete values of 0.8108, 0.8109, and 0.8110. This is because that the Galfit 
software has rounded it final output results to only four decimal places, any difference after 
the four decimal places can not be apparent. 
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Fig. 8.— Optical and Hi properties of the full LSB galaxy sample. Panels respectively 
represent distributions of (a) B-band central surface brightness with a bin of 0.2 magarcsec~‘^, 
(b) B-band absolute magnitude with a bin of 0.2 mag, (c) B-V color with a bin of 0.05 
mag, (d) central surface brightness against apparent magnitude in B band, (e) the common 
logarithm of Hi mass with a bin of 0.1 dex, (f) Hi line velocity width with a bin of 10 kms~^, 
(g)Hl mass against B-band luminosity (measured by our own photometry), and (h) central 
surface brightness against Hi mass . 
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Fig. 9.— Cone diagrams showing the distribution of a.40 Hi sources (blue open circles) and 
the galaxies of our LSB galaxy sample (red open circles) within both the spring and fall sky 
strip covering +24° <decl.<+28°. The diagram shows the volume extending over the full 
ALFALFA velocity range to 18000 km s“^. . 
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Fig. 10.— Distributions of detectable galaxy number counts by SDSS DR7 within a radius 
of 0.5 Mpc (a), 1.0 Mpc (b), and 3.0 Mpc (c) and a redshift gap of z ± 0.04(H-z) for our 
LSB galaxy sample (blue) comparing with all the a.40 galaxies (black). We normalize both 
peaks of blue and black to 1.0 so as to make them easily compared. The grey dashed is 
representing Ngai=8 which is the lower limit of the neighbouring galaxy counts of a galaxy 
candidate in a cluster, and the grey solid line is representing Nga;=16 which is the peak of 
the neighbouring galaxy count distribution for galaxy candidates in clusters. We also show 
distribution of the detectable galaxy number count with a radius of 1.0 Mpc against the Hi 
mass for this LSB galaxy sample in panel (d), in which the grey broken line represents the 
trend between mean Hi mass and mean number count of neighbouring galaxies in 9 Hi mass 
bins with a binsize of 0.5 dex in log M(Hl)/M©. . 
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Fig. 11.— Distribution of the bulge-to-total ratio (B/T) of our LSB galaxy sample. Here, 
the B/T is derived from the parameter “fracDeVr” from the SDSS DR7 catalogue which 
represents the fraction of r-band luminosity contributed by the bulge relative to the lumi¬ 
nosity contributed by all components (bulge -|- disk) of galaxies. The dashed black show the 
bulge-to-total ratios of the 1129 galaxies of our LSB galaxy sample. For a comparison, the 
B/T distribution of the 12423 galaxies of our parent sample is overplotted as solid grey. Both 
histograms are normalized to have a maximum value equal to 1.0 for a clear comparison. 
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Fig. 12.— Galfit fitting by a single component and a combine of two components. The top 
panel shows the exponential profile study for the galaxy AGC 242175 in our LSB galaxy 
sample, and the bottom panel shows the two-component (exponential -|- Sersic) study for it. 
From left to right are respectively the SDSS r-band image, the final model of the galaxy and 
the residual image formed by subtracting the final model from the r-band image. 
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Fig. 13.— Comparison of central disk surface brightnesses of 20 galaxies derived respectively 
from a single component and a combine of two components fitting. The top panel (a) shows 
the comparison of the B-band central disk surface brightnesses derived by a single exponential 
htting, /io(B) (exp) with those derived by a two component of exponential and Sersic (exp 
+ sersic) prohle htting /io(B)(exp+sersic). The bottom panel (b) shows the differences of 
/io(B)(exp+sersic) from /io(B) (exp) versus the bulge-to-total luminosity ratio (B/T). The 
surface brightness is given in units of mag arcsec~‘^. 





- 49 - 

REFERENCES 

Abazajian K. N., Adelman-McCarthy J. K., Agueros M. A., et al. 2009, ApJS, 182, 543 
Adami C. et al. 2006, ApJ, 728, L47 
Adelman-McCarthy K. J. et al. 2006, ApJS, 162, 38 
Adelman-McCarthy K. J. et al. 2008, ApJS, 175, 297 

Beijersbergen M., de Blok W. J. G., van der Hulst J. M. 1999, A&A, 351, 903 
Berlin E. & Arnouts S. 1996, A&AS, 317, 393 

Aihara H., Allende Prieto C., An D., Anderson S. F., Aubourg E. et al. 2011, ApJS, 193, 29 
Blanton M. R. et al. 2005, ApJ, 631, 208 
Bothun G. D. et al. 1997, PASP, 109, 745 
Burkholder V. et al. 2011, AJ, 122, 2318 

Canizares, C. R.,Grindlay, J. E., Hiltner, W. A., Liller, W., McClintock, J. E. 1978, ApJ, 
224, 39 

Gao C., Wu H. 2007, AJ, 133, 1710 
Ceccarelli L. et al. 2012, MNRAS, 426, L6 

Chonis T. S., Martinez-Delgado D., Gabany R. J., Majewski S. R., Hill G. J. et al. 2011, 
AJ, 142, 166 

Courteau S. 1996, ApJS, 103, 363 


Das M. et al. 2009, ApJ, 693, 1300 



de Blok W. J. G., van der Hulst, Bothun G. D. 1995, MNRAS, 274, 235 

de Blok W. J. G., McGaugh S. S., van der Hulst J. M. 1996, MNRAS, 283, 18 

Disney M. & Phillipps S. 1983, MNRAS, 205, 1253 

Duan Z. Y. 2006, AJ, 132, 1581 

Freeman K. G. 1970, ApJ, 160, 811 

Galaz G. et al. 2011, ApJ, 728, 74 

Gao D., Liang Y. G., Liu S. F., Zhong G. H., Glien X. Y. et al. 2010, RAA, 10, 1223 
Geller M. J. et al. 2012, AJ, 143, 102 

Giovanelli R.,Haynes M. P.,Kent B. R., et al. 2005a, AJ, 130, 2598 

Giovanelli R.,Haynes M. P.,Kent B. R., et al. 2005b, AJ, 130, 2613 

Giovanelli R. 2007, Nuovo Gimento B, 122, 1097 

Gunn J. E., et al. 1998, AJ, 116, 3040 

Gunn, J. E., et al. 2006, AJ, 131, 2332 

Haberzettl L. et al. 2007, A&A, 471, 787 

Haynes M. P. 2007, Nuovo Gimento B, 122, 1109 

Haynes M. P., Giovanelli R., et al. 2011, AJ, 142, 170 

He Y. Q., Xia X. Y., Hao G. N., Jing Y. P., Mao S. et al. 2013, ApJ, 773, 37 

Huang S. et al. 2012, AJ, 143, 133 


Hyde B. J., Bernard! M. 2009, MNRAS, 394, 1978 



Impey C. D., Sprayberry D., Irwin M. J., Bothun G. D. 1996, ApJS, 105, 209 
Impey C., Bothun G. 1997, AR&AA, 35, 267 
Impey G. et al. 2001, AJ, 122, 2341 
Infante L 1987, A&A, 183, 177 

Kniazev A.,Grebel E.,Pustilnik S., Pramskij A., Kniazeva T. et al. 2004, AJ, 127, 704 
Kron R. G. 1980, ApJ, 43, 305 

Lauer T. R., Faber S. M., Richstone D., Gebhardt K., Tremaine S. et al. 2007, ApJ, 662, 
808 

Li H. N., Wu H. Gao G., Zhu Y. N. 2007, AJ, 134, 1315 
Liang Y. G., Hu J. Y., Liu F. S., Liu Z. T. 2007, AJ, 134, 759 

Liang Y. G., Zhong G. H., Hammer F., Ghen X. Y., Liu F. S. 2010, MNRAS, 409, 213 
Lin W. P., Zhou X., Burstein D., Windhorst R. A., Ghen J. S. et al. 2003, AJ, 126, 1286 
Lisker T., et al. 2007, ApJ, 660, 1186 

Liu F. S., Xia X. Y., Mao S., Wu H., Deng Z. G. 2008, MNRAS, 385, 23 

Liu Y., Zhou X., Ma J., Wu H., Yang Y. B. et al. 2005, AJ, 129, 2628 

Lupton, R. H., Gunn, J. E., Ivezic, Z., Knapp, G. R., Kent, S., & Yasuda, N. 2001, in ASP 
Gonf. Proc. 238, Astronomical Data Analysis Software and Systems X, ed. F. R. 
Hamden, Jr., F. A. Primini, & H. E. Payne (San Francisco: ASP), 269 

Mao Y. W., Kong X., Lin L. 2014, ApJ, 789,76 


McGaugh S. S. & Bothun G. D. 1994, AJ, 107, 530 



McGaugh S. S., Schombert J. M., Bothun G. D. 1995, AJ, 109, 5 

McGaugh S. S. 1996, MNRAS, 280, 337 

Minchin R. F. et al. 2004, MNRAS, 355, 1303 

Mo H. J. et al. 1998, MNRAS, 295, 319 

Monnier Ragaigne et al. 2003, A&A, 405, 99 

Noguchi M. 2001, ApJ, 555, 289 

O’Neil K. et al. 2000, AJ, 119, 136 

O’Neil K., Bothun G. D., Gornell M. E. 1997a, AJ, 113, 1212 
O’Neil K. et al. 1997b, AJ, 114, 244 

O’Neil K., Bothun G. D., Driel W., Monnier Ragaigne D. 2004, A&A, 428, 823 
Peebles P. et al. 2001, ApJ, 557, 495 

Peng G. Y., Ho L. G., Impey G. D., Rix H. 2002, AJ, 124, 266 

Poggianti B. M., Bridges T. J., Mobasher B., Garter D. et al. 2001, ApJ, 562, 689 

Pustilnik S. A. et al. 2011, MNRAS, 417, 1335 

Roennback J. & Bergvall N. 1995, A&A, 302, 353 

Schlegel D., Finkbeiner D. P., Davis M. 1998, ApJ, 500, 525 

Schombert J. et al. 1992, AJ, 103, 1107 


Schombert J., et al. 2013, AJ, 146, 41 



- 53 - 

Sersic J.L. 1968, Atlas de Galaxias Australes, Observatorio Astronomico, Cordoba, 
Argentina 

Simien, F. & de Vaucouleurs G. 1986, ApJ, 302, 564 

Smith J. A., Tucker D. L., Kent S., Richmond M. W. et al. 2002, AJ, 123, 2121 
Sprayberry D. et al. 1997, AJ, 481, 104 
Stoughton, C., et ah 2002, AJ, 123, 485 
Strauss, M. A., et al. 2002, AJ, 124, 1810 

Trachternach C., Romans D. J., Haberzettl L., Dettmar R. J. 2006, A&A, 458, 341 

van der Hulst J. M. et al. 1993, AJ, 106, 548 

van den Hoek L. B. et al. 2000, A&A, 357, 397 

Vollmer B. et al. 2013, AJ, 145, 36 

Wen Z. L., Han J. L., Liu F. S. 2009, ApJS, 173, 197 

Wu H., Burstein D., Deng Z. G., Zhou X., Shang Z. H., et al. 2002, AJ, 123, 1364 
York D. G., et al. 2000, AJ, 120, 1579 
Zackrisson E. et al. 2005, A&A, 435, 29 

Zheng Z. Y., Shang Z. H., Su H. J., Burstein D., Chen J. S., et al. 1999, AJ, 117, 2757 

Zhong G. H., Liang Y. C., Liu F. S., Hammer F. et al. 2008, MNRAS, 391, 986 

Zhong G. H., Liang Y. C., Hammer F., Chen X. Y., Deng L. C. et al. 2010, A&A, 520, 69 

Zwaan M. A., Briggs F. H., Sprayberry D., Sorar E. 1997, ApJ, 490, 173 



- 54 - 


This manuscript was prepared with the AAS macros v5.2. 



